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Semiclassical Electromagnetic Casimir Self-Energies 

Martin Schaden* 

Rutgers University, 211 Smith Hall, 101 Warren St., Newark, NJ 07102, U.S.A. 

The electromagnetic Casimir energies of a spherical and a cylindrical cavity are analyzed semiclassically. The field 
theoretical self-stress of a spherical cavity with ideal metallic boundary conditions is reproduced to better than 
1%. The subtractions in this case are unambiguous and the good agreement is shows that finite contributions from 
the exterior of the cavity are small. The semiclassical electromagnetic Casimir energy of a cylindrical cavity on the 
other hand vanishes to any order in the real reflection coefficients, whereas the field theoretic Casimir energy of a 
cylindrical boundary that is perfectly metallic and infinitesimally thin is finite and negative|18|. Contrary to the 
spherical case and in agreement with Barton's perturbative analvsis[36|. the subtractions in the spectral density 
for the cylinder are not universal when only the interior modes of are taken into account |4"7|. The Casimir energy 
of a cylindrical cavity therefore depends sensitively on the physical nature of the boundary in the ultraviolet 
whereas the Casimir energy of a spherical one does not. The extension of the semiclassical approach to more 
realistic systems is sketched. 



1. INTRODUCTION 

Demonstrating that the collective interaction 
of atomic systems may have macroscopic con- 
sequences, Casimir obtained the now famous 
attractive force between two neutral metallic 
plates EQ from the boundary conditions they 
impose on the electromagnetic field. Half a 
century later, his prediction has been verified 
experimentally^ to better than 1%. 

Twenty years after Casimir's result for two 
parallel plates, Boycr calculated the zero-point 
energy of an ideal conducting spherical shell[3]. 
Contrary to intuition derived from the attraction 
between two parallel plates, the spherical shell 
tends to expand. Boyer's result has since been 
improved in accuracy and verified by a number 
of field theoretic methods 4l5ltil7l8| - even though 
there may be little hope of observing this effect 
experimentally in the near future[§]. 

Since field theoretic methods require ex- 
plicit or implicit knowledge of cavity frequen- 
cies, they have predominantly been success- 
fully employed to obtain the Casimir energies 
of classically integrable systems. Thus, in ad- 
dition to a spherical cavity, the electromag- 
netic Casimir energies of dielectric slabs |10ll 1 II 2| . 



metallic parallelepipeds 1 3114115116117) and long 
cylinders |7Il8ll9l20l21l!%125| have been com- 
puted in this manner. 

However, most systems are not integrable and 
often cannot even be approximated by such sys- 
tems. It thus is desirable to develop reliable meth- 
ods for estimating the Casimir energies of clas- 
sically non-integrable and even chaotic systems. 
Balian, Bloch and Duplantier calculate Casimir 
energies based on a multiple scattering expan- 
sion to the Green's function |5llMI24| . This ap- 
proach does not require knowledge of the quan- 
tum mechanical spectrum and the geometric ex- 
pansion is exact for sufficiently smooth and ide- 
ally metallic cavities. Ultra-violet divergent con- 
tributions have to be subtracted at every order 
and the relative importance of the finite remain- 
ders is hard to assess a priori. It in practice is 
often difficult to carry the (in principle exact) ex- 
pansion beyond the first few terms without semi- 
classical approximation 2 . In|25| a semiclassical 
method based on Gutzwiller's trace formula 26 
for the response function was proposed to esti- 
mate (finite) Casimir energies. It is appropri- 
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2 Restricting to just two reflections, the Casimir energy 
of an idealized spherical cavity of radius R is of the cor- 
rect sign and quite accurately estimated^ as 3hc/(64R) ~ 
0.0469?ic/i? 



ate for Casimir energies of hyperbolic and chaotic 
systems 25 27 28 with isolated classical periodic 
orbits. The semiclassical approximation to the os- 
cillating part of the spectral density of integrable 
systems was derived in |13l29| and has been used 
to estimate the force between concentric metallic 
cylindrical shellsPOj. The agreement with exact 
results is rather impressive even when the radius 
of the outer cylinder is about four times that of 
the inner one 30 . Although not exact in gen- 
eral, the semiclassical approximation associates 
the finite (Casimir) part of the vacuum energy 
with optical properties of the system. It captures 
aspects of Casimir energies that have been puz- 
zling for some time|31| and appears to be a rea- 
sonable approximation in the experimentally ac- 
cessible regime where these forces are relatively 
large J25I31J . Path integral methods J3H33I34I35| 
in principle allow one to obtain Casimir interac- 
tions between disjoint bodies to arbitrary preci- 
sion. Due to unresolved renormalization issues, 
these methods have so far not been used to study 
the self stress of cavities. The purpose of this ar- 
ticle is to revisit and analyze the Casimir stress 
of (integrable) spherical and cylindrical cavities 
semiclassically and to compare it with field theo- 
retic results. 

The simplicity, transparency and somewhat 
surprising accuracy of the semiclassical approx- 
imation is demonstrated in Boyer's problem J3I4I 
], that is in determining the electromag- 



netic Casimir energy of a spherical cavity with an 
(ideal) metallic boundary. Instead of directly pro- 
ceeding from the semiclassical expressions for the 
oscillating part of the spectral density [T3 29 , our 
starting point will be the dual representation of 
the Casimir energy of integrable systems in terms 
of periodic orbits 26 29 . The subtractions in the 
spectral density that render the Casimir energy 
finite then are apparent and we can better an- 
alyze the semiclassical approximation. Even so, 
the semiclassical analysis of Boyer's problem is an 
order of magnitude simpler than any given pre- 
viously. However, since no bounds are obtained, 
it at present is not possible to judge the accu- 
racy of the result without comparing to field the- 
oretic calculations [B]. It will become rather clear 
though, that the semiclassical analysis is accurate 



enough to infer the sign of the Casimir energy of a 
cavity by geometric arguments if the contribution 
from periodic orbits does not vanish. 

In ref . |3U| periodic orbits in fact were found 
not to contribute to the Casimir energy of a long 
cylindrical cavity. The optical phases that give 
a positive semiclassical Casimir energy for the 
spherical cavity, lead to a vanishing one in the 
cylindrical case. The somewhat intriguing re- 
sult that the Casimir energy of a cylinder |7I2 11241 
36 37 38 vanishes to first order in the reflection 
coefficients thus apparently is readily explained 
by geometric optics. However, the semiclassi- 
cal Casimir energy of a cylindrical cavity van- 
ishes to all orders in the real reflection coefficients 
and thus also vanishes for an ideal metallic cav- 
ity. The discrepancy to the finite field theoretic 
Casimir energy of an idealized, infinitesimally 
thin cylindrical boundary between non-dispersive 
media with the same speed of Iight (18ll9l21j will 
be traced to the presence of a logarithmic diver- 
gence observed by Bartonplo] in his perturbative 
treatment of the non-ideal dilute case. We will 
argue that the exact cancellation of this diver- 
gence in the field theoretic approach is due to the 
infinitesimal thickness of the boundary - interior 
and exterior contributions to the Casimir energy 
in this case depend on only one common scale - 
the radius of the cylinder. 



2. The Dual Picture: Casimir Energies in 
Terms of Periodic Rays 

Integrable systems may be semiclassically 
quantized in terms of periodic paths on invari- 
ant tori|39| - in much the same manner as Bohr 
first quantized the hydrogen atom. Although in 
general not an exact transformation, classical pe- 
riodic orbits on the invariant tori are dual to the 
mode frequencies in the semiclassical sense. Ap- 
plying Poisson's summation formula, the semi- 
classical Casimir energy (SCE) due to a massless 
scalar may be written in terms of classical peri- 
odic orbits |2^9l3T |. 
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The components of the <i-dimensional vector I 
in Eq. QJ are the actions of a set of properly nor- 
malized action-angle variables that describe the 
integrable system. The exponent of the integrand 
in Eq. ||TJ is the classical action (in units of K) of a 
periodic orbit that winds m, times about the i-th 
cycle of the invariant torus. H(T) is the associ- 
ated classical energy and /3 m is the Keller-Maslov 
index (4T)l41| of a class of periodic orbits identi- 
fied by the vector of integers m. The latter is 
a topological quantity that does not depend on 
the actions I. To leading semiclassical order, the 
(primed) sum extends only over those sectors m 
with classical periodic paths of finite action (see 
below). The correspondence in Eq.Q can only 
be argued semiclassically 26 29 and the integrals 
on the RHS therefore should be evaluated in sta- 
tionary phase approximation (sp). 

Contributions to the Casimir energy from high 
frequencies correspond to those from short pe- 
riodic orbits in this dual picture. Divergences 
due to periodic classical paths of vanishing length 
(and thus vanishing action) on the RHS of Eq. JIJ 
are related to ultra-violet divergences of the mode 
sum on the LHS of Eq.QJ. If these divergences 
can be subtracted unambiguously 31 36 42 , the 
dependence of the vacuum energy on macroscopic 
properties of the system is semiclassically repre- 
sented by contributions due to classical periodic 
orbits of finite action only. The primed sum on 
the RHS of Eq.Q indicates this restriction 3 . The 
(divergent) Wcyl contribution to the vacuum en- 
ergy from the m = (0, . . . , 0)-sector in particular 
has to be subtracted. Together with an evalu- 
ation of the integrals in stationary phase, this 
defines the semiclassical Casimir energy (SCE). 
To physically interpret the SCE, one has to con- 
sider the implicit subtractions in the spectral 
density |24I31I36| in greater detail. 



3 This is conceptually not so different from considering only 
the contribution of topologically non-trivial "instanton" 
sectors to the vacuum energy of a field theory. 



3. The Spherical Cavity 

The semiclassical spectrum of a massless scalar 
is exact for a number of manifolds without 
boundary [l3*| and the definition of the SCE by the 
RHS of Eq. QJ coincides with the Casimir energy 
of zeta- function regularization in these cases. It 
also is exact for massless scalar fields satisfying 
periodic-, Neumann- or Dirichlet- boundary con- 
ditions on parallelepipeds 14 16 31 as well as for 
some tessellations of spheres |3 1144145) . In |5S] the 
semiclassical approximation was argued to give 
the leading asymptotic: behavior of the Casimir 
energy whenever the latter diverges as the ratio 
of two relevant lengths vanishes. All these crite- 
ria do not apply to the Casimir self-stress of a 
spherical cavity first considered by Boyerj^J. The 
latter is an integrable system, but the semiclas- 
sical spectrum is only asymptotically correct [T^]. 
There furthermore is no ratio of lengths in which 
one might hope to obtain an asymptotic expan- 
sion. One therefore cannot expect the semiclassi- 
cal approximation to be exact in this case. It nev- 
ertheless turns out to be quite accurate. The SCE 
is obtained by performing the integrals of Eq. (JJJ 
in stationary phase and has a very transparent in- 
terpretation in terms of periodic orbits within the 
cavity only. The sign of the SCE of a spherical 
cavity in particular will be quite trivially estab- 
lished and the good agreement supports the con- 
jecture that the contribution from exterior modes 
mainly serves to cancel the ultra-violet divergence 
of the contribution from interior modes in the 
field theoretic approach 46/E] . The observed dis- 
crepancy of 1% compared to the field-theoretic 
result probably can be attributed to the error in 
the semiclassical estimate of low-lying eigenval- 
ues of the Laplace operator - which is of similar 
size. Since boundary conditions for the electro- 
magnetic field never are ideal, many corrections 
of similar and even greater magnitude would be 
required in any realistic situation. 

The electromagnetic Casimir energy of any 
closed cavity with a smooth and perfectly metal- 
lic boundary may be decomposed into the con- 
tribution from two massless scalar fields - one 
satisfying Dirichlet's, the other satisfying Neu- 
mann's boundary condition on the surface |31). 



Because the surface is ideally metallic, contri- 
butions to the spectral density from arbitrarily 
short closed paths that reflect off (either side of) 
the smooth surface cancel each other. Semiclassi- 
cally there is no (potentially divergent) local con- 
tribution to the Casimir energy from such an ide- 
alized surface in the electromagnetic case - the 
local surface tension in fact vanishes for an in- 
finitesimally thin surface |24!48| . Note that this 
cancellation is quite special for the electromag- 
netic field and metallic boundaries. It in general 
does not occur for a massless scalar field satisfying 
Neumann or Dirichlct boundary conditions on a 
spherical surface of arbitrary thickness or in even 
dimensions 46 49 5(JI51| . However, it was shown 
in |3B] that the Casimir energy due to a scalar 
is finite for infinitesimally thin spherical cavities 
in odd dimensional spaces. The following argu- 
ment illustrates the absence of ultraviolet diver- 
gent contributions proportional to the "area" of 
any infinitesimally thin even-dimensional surface 
if the speed of light in the surrounding medium 
is constant. Barring other scales and with a 
dimensionless cutoff 4 , the local contribution to 
the surface divergence from a small ((d — 1)- 
dimensional) surface element dA is proportional 
to hcdAY,ifi{{Rj/Rk})/Ri, where Ri are (lo- 
cal) radii of curvature and the /, are dimension- 
less functions of their ratios only. The surface 
divergence semiclassically arises due to arbitrary 
short closed optical paths on either side of the in- 
finitesimally thin boundary that touch the surface 
at a single point only. But the radii of curvature 
on either side of a point on the infinitesimally 
thin surface are of equal magnitude and of oppo- 
site sign. This implies that divergent contribu- 
tions from the two arbitrarily short closed optical 
paths with a common point on the surface that 
lie on either of its sides cancel precisely for odd d 7 
whereas they are equal in even d - provided the 
boundary conditions and speed of light are the 
same for both sides of the surface. One therefore 
expects the scalar Casimir energy of such cavi- 



4 In dimensional regularization this cutoff is just the 
(small) deviation from integer d and in the semiclassical 
context it is the ratio of the cutoff in the length of the 
shortest optical paths to some characteristic (fixed) length 
in the problem. 



ties to be analytic in the vicinity of odd d, but 
to in general diverge for even dimensional spaces. 
The finite Casimir energy of thin spherical cavi- 
ties in d dimensions 0H1 as wei l as the perturba- 
tive analysis of such surface divergences |SJ agree 
with this semiclassical analysis. In three dimen- 
sions, the Casimir energy due to a scalar field 
was found to be finite |19l47l52| also for an in- 
finitesimally thin cylindrical cavity. The above 
argument is quite general and indicates that sur- 
face divergences in will cancel locally only for in- 
finitesimally thin (and sufficiently smooth) even- 
dimensional boundaries when the speed of light is 
the same on either side. This cancellation of sur- 
face divergences is due to the infinitesimal thin- 
ness of the boundary and does not depend on 
the precise nature of the (frequency-independent) 
boundary condition. [It is quite different from the 
semiclassical cancellation of surface divergences of 
the electromagnetic field on ideal metallic bound- 
aries in three dimensions. The latter occurs due 
to the different boundary conditions for the two 
scalar modes and does not depend on the thick- 
ness of the metallic surface.] 

The only subtraction in the spectral density re- 
quired for a finite Casimir energy in the electro- 
magnetic case with idealized metallic boundary 
conditions thus is the Weyl contribution propor- 
tional to the volume of the sphere. The latter 
corresponds to ignoring the m = (0, 0, 0) con- 
tribution to the sum in Eq.JQ). The remain- 
ing difficulty in calculating the SCE of an inte- 
grable system is a convenient choice of action- 
angle variables. For a massless scalar in three 
dimensions satisfying boundary conditions with 
spherical symmetry, an obvious set of actions is 
the magnitude of angular momentum, I2 = L, one 
of the components of angular momentum I3 = L z 
and an action I\ associated with the radial degree 
of freedom. 

Since the azimuthal angle of any classical orbit 
is constant, the energy E = H{I\,l2) of a mass- 
less particle in a spherical cavity of radius R does 
not depend on I3 = L z . In terms of the previous 
choice of actions, the classical energy is implicitly 



given by, 
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The branches of the square root and inverse cosine 
in Eq.JSJ) are chosen so that I\ is positive. It is 
convenient to introduce dimcnsionless variables 



X = 2ER/(hc) and z = d 2 /(ER), 



(3) 



for the total energy (in units of Tic/(2R)) and 
the angular momentum (in units of ER/c) of an 
orbit. Note that z 6 [0, 1] and that the semi- 
classical regime formally corresponds to A ^> 1, 
i.e. to wavelengths that are much shorter than 
the dimensions of the cavity. With the help 
of Eq.J2J) and the definitions of Eq.lQ, the semi- 
classical expression in Eq.QJ for the Casimir en- 
ergy of a massless scalar field satisfying Neumann 
or Dirichlet boundary conditions on a spherical 
surface becomes, 
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The integral over I3 has here been performed 
in stationary phase approximation. Because the 
Hamiltonian does not depend on I3 , only periodic 
orbits with 771.3 = contribute. Since — 12 < ^3 < 
I2, one has that f (H3 — 2I 2 = fiXz. The fac- 
tor 2/2 accounts for the 2(1 + l/2)-degeneracy of 
states with angular momentum L = h(l + 1/2) = 
I-2\. By taking (4 times) the real part in Eq.Q 
one can restrict the summations to non-negative 
integers and choose the positive branch of the 
square root- and inverse cosine- functions in the 
exponent 5 . The Keller-Maslov index f3(n, m) of a 
classical sector depends on whether Neumann or 
Dirichlet boundary conditions are satisfied on the 
spherical shell and will shortly be determined. 

For positive integers 777, and 77, the phase of the 
integrand in Eq.Q is stationary at z = z(n, 777) e 

5 The primed sum now implies half the summand if one of 
the integers vanishes as well as the absence of the m = 
n = term. 



[0, 1] where, 

= —77 arccos(z) + 777,7r 

=>■ £(77,777) = cos(rmr/n), n > 2t77 > 1 . (5) 

Restrictions on the values of 777 and 77 arise be- 
cause arccos(z) <E [0, 7r/2] on the chosen branch. 
The phase is stationary at classically allowed 
points only for sectors with 77 > 2m > 1. Semi- 
classical contributions to the integrals of other 
sectors arise due to the endpoints of the z- 
integration at z = and z = 1 only. Such 
" diffractive" contributions are of sub-leading or- 
der in the asymptotic expansion of the spectral 
density for large A. Note that 777 —* m+n amounts 
to the choice of another branch of the inverse co- 
sine. 

The classical action in sectors with stationary 
points is, 

S c i(n,m) — hXn sm(mir/n) (6) 

= (E/c)2nRsm(rmr/n) = (E/c)L(n,m) , 

where L(n, m) is the total length of the classical 
orbit. Some of these periodic orbits are shown in 
Fig. 1. The integer 777 in Eq.@ gives the number 
of times an orbit circles the origin. The integer 
77 > 1 gives the number of times an orbit touches 
the spherical surface. As indicated in Fig. 1, the 
set of classical periodic orbits in the (77, m)-sector 
form a caustic surface and a double covering is 
required for a unique phase-space description [4()|. 
The two sheets are joined at the inner caustic 
[indicated by a dashed circle in Fig. 1] and at 
the outer spherical shell of radius R. Every orbit 
that passes the spherical shell 77 times also passes 
through the caustic 77 times. The cross-section 
of a bundle of rays is reduced to a point at the 
spherical caustic surface. The caustic thus is of 
second order and is associated with a phase loss of 
7T every time it is crossed. At each specular reflec- 
tion off the outer shell, Dirichlet boundary condi- 
tions require an additional phase loss of 7r whereas 
there is no change in phase for Neumann bound- 
ary conditions. Altogether the Keller-Maslov in- 
dex of sector (77, 777) depends on 77 only and is 
given by, 



j3(n, m) = 



0, for Dirichlet b.c. 
2?7, for Neumann b.c. 



(7) 



Since the electromagnetic Casimir energy of a 
cavity with a smooth surface on which (ideal) 
metallic boundary conditions hold can be viewed 
as due to two massless scalar fields, one satisfy- 
ing Dirichlet and the other Neumann boundary 
conditions |24| on this surface, only sectors (n, to) 
with even n — 2k > 2to > 2 contribute j24j to the 
SCE in leading order of the asymptotic expansion 
for large A. 




Fig.l: Classical periodic rays within a ball or a solid 
cylinder, a) The shortest primitive rays with winding 
numbers (n,m) g {(2, 1), (3, 1), (4, 1)}. b) Primitive 
rays to winding numbers (n,m) — (5,1) and (5,2). 
Caustic surfaces are shown as thin circles. The dashed 
part of any trajectory is on one sheet and its solid 
part on the other of a two-sheeted covering space. 
The "phase space" of the (5, 2) sector is indicated by 
the hatched area. Note that caustics are of 2 nd order 
for a spherical cavity but of 1 st order for a cylindrical 



Note that sectors with m = or n = have 
vanishing classical action and do not contribute 
to the SCE. Eq.(0| implies that extremal paths 
in the (n > 0, to = 0) sectors have maximal angu- 
lar momentum 1 = z = lc/(ER). These are great 
circles that are wholly within the spherical shell 
in a plane perpendicular to the one under con- 
sideration, i.e. classical orbits with I3 = L z = 0. 
Because the measure of the z-integral vanishes at 
z — 1 like y/1 — z these classical paths are ex- 
tremal but not stationary. This can also be seen 
by expanding the exponent in Eq.Q about the 
stationary point z(n,m). For (n > 0,to > 0) the 
second derivative of the classical action is finite 



at z(n, m), 

3 Ti 
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,(8) 



whereas it diverges in sectors with to = 0. The 
behavior of the exponent for z ~ 1 in this case is, 



vT^I2--zarccos(z) = ^fl(i-zf' 2 +0({i-z 



,5/2^ 



.(9) 



Quadratic fluctuations about the classical orbit 
with to = thus have vanishing width and these 
sectors do not contribute in stationary phase ap- 
proximation. To leading semiclassical accuracy, 
the Casimir energy of a spherical cavity with an 
ideal metallic boundary therefore is. 
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This semiclassical estimate is only about 
1% larger than the best numerical value 
0. 04617. ..hc/R for the electromagnetic Casimir 
energy of a spherical cavity with an infinitesimally 
thin metallic surface. The error is of the same 
order as that of limiting the multiple reflection 
expansion to just two reflections, which gives[5] 
£bj^ ~ 0.0469?ic/i?. Note that the contribution 
from the (2k, k) sectors had to be considered sep- 
arately in Eq. I|10l) since the measure dzz vanishes 
at the stationary point z(2k, k) — cos(7r/2) = 
of the integrand, which is an endpoint of the in- 
tegration domain. As can be seen in Fig. la), 
the classical rays of (2k, fc)-sectors go back and 
forth between antipodes of the cavity and pass 
through its center - they have angular momen- 
tum L = = z(2k,k). 

The shortest primitive orbits give somewhat 
less than half (1/(16tt) - 0.02) of the total SCE 
of the spherical cavity - much less than the 92% 
they contribute to the Casimir energy of paral- 
lel plates. The main reason is that contributions 



only drop off as 1/fc 2 rather than like 1/fc 4 as 
for parallel plates. The length of an orbit in 
the (4, l)-sector (the inscribed square in Fig. la) 
furthermore is just a factor of v2 longer than 
a (2, l)-orbit [which in turn is a factor of 1/V2 
shorter than a (4, 2)-orbit]. To estimate the mag- 
nitude of the contribution from any particular 
sector one has to take the available phase space 
as well as the ray's length into account. Thus, al- 
though the length of a (2k, l)-orbit tends to 2irR 
for k — > oo, the associated phase-space (essen- 
tially given by the volume of the shell between 
the boundary of the cavity and the inner caustic) 
decreases like 1/fc 2 . This accounts for the rela- 
tively slow convergence of the sum in Eq. i|l(J|) . To 
achieve an accuracy of 10 -5 , the first 50 terms of 
the sum were evaluated explicitly and the remain- 
ing contribution was estimated using Richard- 
son's extrapolation method. 

4. The Cylindrical Cavity 

The example of a spherical cavity shows that 
the SCE in some instances is surprisingly ac- 
curate. However, there evidently are systems 
without periodic classical orbits, such as the two 
perpendicular planes investigated in 33] , or the 
Casimir pendulum of|53). None of these systems 
is integrable, and although there are no stationary 
periodic classical rays, periodic rays of extremal 
(shortest) length do exist. Semiclassically, such 
extremal periodic rays are associated with diffrac- 
tion 54 55 . The inclusion of diffractive contribu- 
tions in the semiclassical estimate of Casimir en- 
ergies has so far only been attempted for a system 
of spheres 56 . It will become evident below that 
diffractive contributions also play a central role 
in the Casimir energy of a cylindrical cavity. 

The Casimir energy of a dilute cylindrical gas 
of atoms was found to vanish in|57). A num- 
ber of calculations have confirmed that there is 
no contribution up to second order in the reflec- 
tion coefficients for dielectrics 36 37 38 and for 
media where the speed of light on either side of 
an infinitesimally thin cylindrical boundary is the 
same j7l2()l21l22l24) . Balian and Duplantier even 
conjectured that the Casimir energy of an ideal 
metallic cylindrical cavity may vanish 24 to all 



orders of the multiple reflection expansion. The 
non-vanishing Casimir energy of an ideal metal- 
lic cylindrical cavity [T^j was reanalyzed in the 
framework of zeta-function regularization. It was 
confirmed that the Casimir energy of an ideal 
metallic cylinder only vanishes to leading order 
and that higher orders in the reflection coeffi- 
cients all give a non- vanishing contribution^^]. 
However, some mathematical prowess is required 
to analytically prove the lowest order cancella- 
tion in the field-theoretic approach 37 38 . That 
a number of separate contributions should con- 
spire to a null result without apparent physical 
reason has been considered by many as somewhat 
"mysterious" |7ll2| . The suspicion that this can- 
cellation could be a purely geometrical effect, is 
nourished by the fact that the finite part of the 
pair-wise Van DerWaals interaction energy of a 
dilute gas of atoms vanishes for a cylinder |36l57j 
but not for other geometries. In 30 it was found 
that the semiclassical contribution due to peri- 
odic rays also vanishes for a cylinder. However, a 
careful perturbative analysis reveals that the in- 
teraction energy of any real dilute cylindrical gas 
of atoms includes a logarithmic divergence in ad- 
dition to divergent contributions proportional to 
the volume and surface area of the cylinder |36|. 
The subtraction of this logarithmic divergence 
generally is ambiguous and the Casimir energy of 
a cylindrical cavity depends sensitively on prop- 
erties of its boundary 58 in the ultraviolet. A 
particular boundary (say an infinitesimally thin 
cylindrical shell separating media with the same 
speed of light) thus may have a finite (negative) 
Casimir energy, whereas a very small modifica- 
tion of this boundary (say finite thickness) leads 
to a logarithmic divergence. 

The calculation below supports this possibility. 
It is already known fromPO] that the semiclassi- 
cal contribution to the Casimir energy vanishes 
for a cylinder. The SCE of a cylindrical cavity 
in fact vanishes to all orders in the reflection co- 
efficients for the same reason that the SCE of a 
spherical cavity is positive - due to relatively ob- 
vious optical phases. The semiclassical point of 
view thus gives a straightforward and physically 
acceptable explanation for otherwise mysterious 
cancellations. It also indicates that any addi- 



tional phase change at the boundary will destroy 
this delicate mechanism. The finite electromag- 
netic Casimir energy[l8 of a cylinder with ide- 
alized metallic boundary conditions on the other 
hand is more difficult to explain semiclassically. 
However, contrary to a spherical cavity and in 
agreement with the perturbative result of 36 , the 
semiclassical expression for the Casimir energy of 
a cylindrical cavity also is logarithmic divergent. 
This is due to " diffractive" end-point contribu- 
tions that are ignored in stationary phase approx- 
imation. There is reason to believe fl7!52| that the 
subtraction of this logarithmic divergence by the 
contribution from "exterior" modes gives a finite 
Casimir energy for an ideal metallic and infinites- 
imally thin cylinder 18 . 

To exhibit these effects we revisit the calcula- 
tion of the electromagnetic SCE of a long cylin- 
drical cavity, or rather of a very thin torus with 
one perimeter L that is much larger than the 
other, L ^> 2irR. This is an integrable system. 
In the limit R/L — ► 0, the only classical trajec- 
tories of relevance are again those of Fig. 1 and 
the SCE of a long cylindrical cavity can be ob- 
tained along similar lines as that of a spherical 
one - with some important modifications. Due to 
the toroidal symmetry of the (long) cylinder, the 
third action I3 = Lp^/ (2tt) in this case is propor- 
tional to the conserved momentum pl along the 
axis of the (thin) cylinder and in Eq. (J2J the en- 
ergy E must be replaced by \/E 2 — i^-nd^/V) 2 . 
The second action furthermore is the angular mo- 
mentum rather than just its magnitude. It again 
is convenient to consider dimensionless quantities 
for the fraction — 1 < x < 1 of the total momen- 
tum along the axis of the cylinder, for the ratio 
— 1 < z < 1 of the angular momentum to the 
maximal possible angular momentum of a photon 
within the cavity and for its energy < A < 00 
in units of hc/(2R), 



A = 2ER/{Tic) 



ch/{ER) 
VT^ 2 • 



2ttcI 3 
EL 



(11) 



Proceeding as in the spherical case, the semiclas- 
sical expression in Eq.Q for the SCE of a mass- 
less scalar field satisfying Neumann or Dirichlct 
boundary conditions on a cylindrical surface be- 



comes, 



£cvl — 



-'cyl 



hcL s—^ 

16ir 2 R 2 ^ 

m,n>0 
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(12) 

The contribution from periodic orbits that wind 
around the perimeter of the torus is negligible 
in the R/L — > limit and has been omitted 
in Eq. (jT2")l . The phase of the integrand in Eq. lfT2T) 
is stationary at x = (corresponding to pl = 0) 
and z(n,m) given in Eq.Q. Since the domain 
of integration for the z-variable differs from the 
spherical case, sectors with 1 < m < n — 1 have 
non-trivial stationary points. The classical action 
of an (n, m)-sector is the same as for the spheri- 
cal cavity and is given by Eq. © . The fluctuations 
about such a classical ray on the other hand are 
quite different for cylindrical and spherical cavi- 
ties. To quadratic order in the fluctuations about 
the stationary point x = 0, z(n, m), the action for 
the cylinder is 



On 



x 2 . mir (z — z(n,m)) 2 
(1 -T )Sm - + 2siU 



(13) 



The unconstrained Gaussian integrals over z — 
z(n,m) and x result in a factor of 27r/(nA) in 
stationary phase approximation. Note that the 
phases of ±7r/4 associated with the two Gaussian 
integrals cancel in this case. Performing also the 
integral over A in Ea. l|12l) finally gives, 



^cyl 



hcL 
AnR 2 



EE^ : 



-i^(3{n,m) 



n 4 sm 2 « 



(14) 



The crucial difference to a spherical cavity is the 
phase factor of —i. It arises because the fluc- 
tuations of a cylindrical system have one fewer 
zero-mode than those of a spherical one 6 . The 



6 The Hamiltonian of a spherical cavity does not depend on 
I3 oc L z , whereas it does depend on ^3 oc pl for the cylin- 
drical cavity. The 2x2 Hessian matrix Hij = d 2 H/dIidIj 
with 3 > i, j > 1 has one zero mode for a spherical cavity, 
but none for the cylindrical geometry. This difference in 
zero modes implies |3"T1 an additional phase loss of n/2 for 
the periodic rays of a cylindrical cavity. 



additional phase loss of 7r/2 ultimately is respon- 
sible for the vanishing of the SCE of a cylindri- 
cal cavity. To verify this we only need to com- 
pute the Keller-Maslov index f3(n, m) for Neu- 
mann and Dirichlet boundary conditions. The 
caustics of the cylindrical cavity are of first order 
rather than second: the cross-section of a bundle 
of rays becomes one-dimensional at the caustic - 
it is focussed to a line rather than a point. Taking 
into account the phase retardation by 7r/2 every 
time a ray passes a first order caustic, the analo- 
gous result to Eq.Q for a cylindrical cavity is, 



f3(n,m) 



3n, for Dirichlet b.c. 
n, for Neumann b.c. 



(15) 



Contributions from paths with Neumann and 
Dirichlet boundary conditions and an odd num- 
ber of reflections cancel each other and, as for the 
spherical cavity, only sectors to even n = 2k = 
2,4,... contribute to the electromagnetic SCE 
[for smooth metallic cavities this is quite gener- 
ally so|23]]. Summing contributions to the elec- 
tromagnetic Casimir energy from the two scalars 
in Eci. (|14J) then gives the null result 
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(16) 



In Eq. (|16[) every periodic orbit gives a vanish- 
ing contribution to the SCE of a cylindrical cav- 
ity. The cancellation evidently depends on a del- 
icate relation between the optical phases. It is 
interesting that a small additional phase loss at 
each reflection off the surface results in a neg- 
ative SCE for a cylindrical cavity, but that the 
Casimir energy vanishes as long as the above 
phase relations hold - even if the magnitude 
of the reflection coefficients is less than unity. 
The SCE in this sense is in line with previous 
results jSEEZEHEIj for the Casimir energy of a 
dilute dielectric cylinder, and in fact supports 
the conjecture of Balian and Duplantier in|24). 
The non- vanishing Casimir energy of a cylindrical 
cavity with an infinitesimally thin ideal metallic 
boundary on the other hand is not so easily ex- 
plained by this semiclassical point of view. 

Some insight is gained by noting that the con- 
tribution of any sector to the SCE of a cylin- 



drical cavity in Ea. i|12f) - even sectors with non- 
trivial periodic classical paths - diverges. This is 
in marked contrast to the spherical case, where 
the contribution from sectors with non-trivial pe- 
riodic classical paths (characterized by n > 2m > 
1) is finite. The divergence is readily made ex- 
plicit by scaling X^/l — x 2 — > A in the integral 
of Ea.(|12p. Without ultraviolet cutoff, the result- 
ing x-integral in this case formally gives the fac- 
tor, 



dx 



-i(l-* 2 ) 2 



oo 



(17) 



which diverges due to the behavior of the in- 
tegrand near the endpoints of the integral at 
x = ±1. It may be regulated by introducing an 
ultraviolet cutoff 51 of some sort for the energy 
integral [that is in the integral over A]. As may 
be seen from Eq. (|17|l . the regulated integral will 
always include terms that are logarithmically di- 
vergent as n — > oo. The subtraction of a logarith- 
mic divergence depends on details of the cutoff 
and thus is sensitive to ultraviolet properties of 
the boundary j42j. The evaluation of (in principle 
divergent) integrals in stationary phase could be 
considered one way of subtracting this divergence. 
Because the divergence is logarithmic, the sub- 
traction is by no means unique in this case. The 
presence of a logarithmic divergence for cylindri- 
cal cavities was first emphasized by Barton|5B] 
in his perturbative treatment of a dilute gas of 
atoms. It also is evident in the contribution from 
interior modes to the Casimir energy of an ideally 
metallic cylinder 47 . 

The foregoing is compatible with previous 
results [T8l2H7| that the Casimir energy of a cylin- 
drical cavity is finite if the speed of light inside 
and outside its infinitesimally thin boundary are 
the same. It for instance is negative for idealized 
metallic boundary conditions ^Hj- The Casimir 
energy in this case apparently does not suffer from 
any logarithmic divergences (or equivalently, from 
any pole ambiguities in zeta function regulariza- 
tion). The Casimir energy is finite for the in- 
finitesimally thin boundary, because the logarith- 
mic divergent contribution from interior modes 
is precisely cancelled by the similarly logarith- 
mic divergent contribution from exterior modes. 
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Since the boundary is infinitesimally thin and the 
speed of light on both sides of the boundary is 
the the same, a precise cancellation is possible. 
The divergence reappears for a dielectric cavity 
in vacuum with a lower speed of light in the 
dielectric [SSj. This occurs for a spherical as well 
as for a cylindrical cavity, but with an impor- 
tant difference: the divergence in the spherical 
case is not logarithmic and may be unambigu- 
ously subtracted [5]. The subtraction of the log- 
arithmic divergence in the Casimir energy of the 
cylindrical cavity on the other hand requires the 
introduction of some energy scale that describes 
ultraviolet properties of the boundary. An anal- 
ogous problem would be encountered for an ideal 
metallic boundary of finite thickness 07] and in 
fact for almost any small deviation from an ide- 
alized and infinitesimally thin cylindrical bound- 
ary between two media with identical speed of 
light. Paradoxically, defining the Casimir energy 
of an cylindrical cavity in a manner that does not 
depend on detailed ultraviolet properties of its 
boundary appears all but impossible. 

It perhaps is worth mentioning in this regard 
that the Casimir energy of a massless scalar 
excitation on the two-dimensional spherical or 
toroidal boundaries is well-defined. For a spher- 
ical shell and a very thin torus, this Casimir en- 
ergy has the same dependence on the dimensions 
as the Casimir energies of the corresponding cav- 
ities. For a two-sphere (52) and a very thin torus 
T2 with L ^S> 2ttR these Casimir energies are 



e 8a = 

£t 2 = - 



(18) 



HcL 



C(3) 



-0.0097. 



HcL 



Att 3 R 2 

Note that these Casimir energies of a mass- 
less scalar on two-dimensional spherical and 
toroidal surfaces are exactly reproduced 
semiclassicallv pril4.S| . The presence of scalar 
surface modes therefore does not change the 
Casimir energy of a spherical cavity but could 
very well contribute to that of a cylindrical one. 
The Casimir energy of a massless degree of free- 
dom on a torus not only is of the same form, 
but also of the same sign and order of magnitude 
as the Casimir energy of an ideal metallic cylin- 
drical cavity [TBI 19j . Such a contribution from 



massless surface modes thus might be important 
for a cylindrical cavity and would be difficult 
to separate from the contribution due to cavity 
modes. 



5. Discussion 

The semiclassical approximation to the Casimir 
energy of a cavity to leading order includes only 
contributions from quadratic fluctuations about 
stationary periodic classical rays. Since all peri- 
odic rays lie in the interior, the SCE of a con- 
cave cavity to leading order depends on the exte- 
rior only indirectly through reflection coefficients. 
Periodic classical rays furthermore are of finite 
length. Their contribution to the Casimir energy 
thus is ultraviolet finite. However, this approx- 
imation is sensible only if UV-divergent contri- 
butions to the vacuum energy can be subtracted 
unambiguously from the spectral density. Loga- 
rithmically divergent contributions to the vacuum 
energy require a subtraction scale 05] . The latter 
is a clear indication that the subtraction cannot 
be universal since it depends sensitively on the 
UV-properties of the boundary. Small changes 
in the boundary conditions in this case do not 
necessarily correspond to small changes in the 
Casimir energy. The local properties of a bound- 
ary apparently include its thickness: whereas the 
Casimir energy of a cylindrical cavity with an 
ideal and infinitesimally thin metallic boundary 
is finite ^Hl to any order in the (real) reflection 
coefficients \Hj\ , a logarithmic dependence on the 
cutoff appears in more realistic situations J36I47] . 
An ambiguous subtraction also is required in the 
semiclassical approximation. The absence of any 
logarithmic divergence for the infinitesimally thin 
boundary apparently is due to a cancellation by 
exterior modes. Such a cancellation of logarith- 
mic singularities can occur when exterior and in- 
terior modes depend on precisely the same scale, 
the radius R of the cylindrical cavity in this case. 
Although the two logarithmic divergences (each 
proportional to hcL/R 2 for dimensional reasons) 
cancel in the idealized situation, they would not 
do so for a boundary of finite thickness. 

We considered only the semiclassical Casimir 
energy (SCE) of a spherical and of a toroidal cav- 
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ity with ideal metallic boundary conditions, that 
is, with real reflection coefficients of unit magni- 
tude. These are integrable systems and the SCE 
was obtained from the "dual" description of the 
spectral density in terms of periodic paths on in- 
variant tori 26 29 . The winding numbers of a 
periodic orbit are dual to the quantum numbers 
of a mode. In stationary phase approximation 
the SCE of a spherical cavity is positive and re- 
produces the field theoretic value for an infinites- 
imally thin metallic boundary to within 1%. The 
calculation is rather short and straightforward 
and leads to the convergent sum of Eq. I jl()(l . Each 
term in this sum may be interpreted as the con- 
tribution from a class of periodic rays. A few of 
the shorter primitive periodic rays are depicted in 
Fig. 1. The contribution from any sector is finite 
in this case. 

The contribution from periodic orbits to the 
SCE of a cylindrical cavity with an ideal metallic 
boundary on the other hand vanishes to all orders 
in the number of reflections 30 . This occurs due 
to an overall phase change by an odd multiple 
of ir/2 for any classical periodic ray. Restrict- 
ing to just two reflections, this null result agrees 
with field theoretic calculations for infinitesimally 
thin metallic boundaries |7I21I23I24I5D] . The van- 
ishing SCE appears to support the conjecture of 
Balian and Duplantier that the Casimir energy 
of a metallic cylindrical cavity may vanish. How- 
ever, contrary to the spherical case, the contribu- 
tions of any classical sector to the SCE of a cylin- 
drical cavity diverges. Without subtraction of the 
UV-divergent part, the (finite) semi-classical con- 
tribution to the vacuum energy we obtained is not 
very meaningful. Unfortunately the divergence of 
the integral in Eq. i|17|) includes a logarithmic de- 
pendence on the cutoff. The subtraction of UV- 
divergent contributions to the Casimir energy of 
a cylinder thus is sensitive to a scale and can- 
not be achieved in a universal fashion. The log- 
arithmic dependence on the cutoff was first ob- 
served by Barton[3B] in his perturbative calcu- 
lation of the vacuum energy for a (dilute) gas 
of cylindrical shape to lowest order in the fine 
structure constant. Semiclassically this would 
also correspond to considering the contribution 
from rays with only two reflections (n = 2). That 



the UV-subtractions are fragile and depend cru- 
cially on the UV-properties of the boundary is 
also observed when the speed of light within and 
outside an infinitesimally thin cylindrical bound- 
ary differ [5S]. In the electromagnetic case, the 
logarithmic divergences of exterior and interior 
contributions to the vacuum energy of a metallic 
cylinder cancel for an infinitesimally thin metallic 
boundary 47 . However, they in general cannot 
be unambiguously subtracted 52 . 

These examples of a spherical- and cylindrical 
cavity show that the SCE is quite reasonable and 
is rather simple to calculate when the Casimir 
energy is robust, that is, when the necessary sub- 
tractions depend on global properties of the sys- 
tem only and do not require the introduction of 
an additional scale. This apparently is not possi- 
ble if the divergence is logarithmic |42l3l] . The 
classical periodic paths that contribute to the 
SCE of a concave cavity in stationary phase ap- 
proximation lie entirely within the cavity. Their 
contribution depends on the exterior of the cav- 
ity through reflection coefficients only. It has 
been argued for some time that a Casimir en- 
ergy obtained without explicit inclusion of exte- 
rior modes (as for a parallelepiped 14 16 l7j) is 
all but meaningless 221- The criterion favored 
here 36 31 considers any definition of a Casimir 
energy reasonable (and in principle physically re- 
alizable) in which the UV-divergences of the vac- 
uum energy have been subtracted without refer- 
ence to local properties of the boundary. The 
subtraction may (and in general will) include di- 
vergent contributions from exterior modes. The 
Casimir energy of a parallelepiped can be consid- 
ered a case in point: as Power |6D] did for just two 
slabs, one can always assemble parallelepipeds to 
a cube with fixed dimensions - the Casimir en- 
ergy of an individual Darallelepiped |14ll6ll7j in 
this case reflects changes in the vacuum energy of 
the whole cube as the three dividing planes are 
moved adiabatically. By moving interior surfaces 
of the cube (that in principle could have finite 
thickness), one measures only the finite part of its 
vacuum energy that depends on the dimensions 
of the individual parallelepipeds. By contrast, 
it is difficult to imagine that global changes in 
a vacuum energy are measurable (or even physi- 



12 



cally relevant) if their finiteness depends critically 
on local characteristics of the system [3TJ. Per- 
haps somewhat surprisingly, the electromagnetic 
Casimir energy due to the interior of a very long 
cylindrical cavity does not appear to be robust in 
this sense, whereas the electromagnetic Casimir 
energy due to the interior of a spherical cavity is. 
Apart from relating Casimir energies to opti- 
cal properties, one of the advantages of a semi- 
classical description would be the possibility to 
model more realistic (but robust) physical sys- 
tems. The previous considerations are readily ex- 
tended to dielectrics by using appropriate com- 
plex and in general frequency-dependent reflec- 
tion coefficients. In the case of dielectric slabs 
Milton has shown [H| that Lifshitz's theory [TOUT] 
may be reproduced in this manner. Finite tem- 
perature is incorporated |6 l| by allowing periodic 
rays to also wrap around a fictitious periodic ex- 
tra dimension of circumference hc/{kT). Finite 
temperature corrections thus are small if some 
classical periodic paths are much shorter than this 
circumference. At room temperature the length 
of a periodic ray increases by about 7.6 microns 
every time it winds about the temperature direc- 
tion. Temperature corrections therefore are tiny 
for most nanometer scale experiments 7 but could 
be of greater interest in some astrophysical con- 
siderations (3°K = 1mm). Corrections due to 
surface roughness generally will be more impor- 
tant in technological applications. Many classical 
models for diffuse reflection from rough surfaces 
exist and Lambert's Law is easily incorporated in 
the semiclassical approach by appropriate reflec- 
tion coefficients. The dependence on the wave- 
length perhaps can be modelled by a term of the 
action that accounts for stochastic fluctuations in 
the length of a classical periodic orbit upon reflec- 
tion from rough surfaces. Apart from an average 
change in length, this leads to a damping term 
of the form — (ALE/hc) 2 /2 in the classical ac- 
tion, where (AL) 2 is the variance in the length of 
the periodic orbit. Assuming that this variance 
is itself proportional to the length of the orbit, 
surface roughness can semiclassically perhaps be 



modelled by the modified dispersion 
cp(E) =E + ieE 2 /(hc) , 



(19) 



7 The claim that this correction has been measured to 
sufficient accuracy|21 to distinguish between different ap- 
proaches has recently been disputed 62 . 



where e is a typical length scale for the (stochas- 
tic) roughness of the surface. The predominant 
effect of the modified dispersion of Eq. IJ19J1 is that 
contributions to the Casimir energy from wave 
lengths A <C £ are very much suppressed. A sim- 
ilar conclusion may be drawn from a recent and 
considerably more sophisticated analysis [05]. 

Acknowledgements: I am indebted to Larry 
Spruch for advice and numerous discussions and 
would like to thank K. A. Milton for a critical 
reading of the manuscript. This work was par- 
tially supported by the National Science Founda- 
tion under Grant No. PHY0555580. 

REFERENCES 

1. H. B. G. Casimir, Kon. Ned. Akad. Wetensch. 
Proc. 51, 793 (1948). 

2. S. K. Lamoreaux, Phys. Rev. Lett. 78, 5 
(1997); U. Mohideen and A. Roy, Phys. 
Rev. Lett. 81, 4549 (1998); T. Ederth, 
Phys. Rev. A 62, 062104 (2000); G. Bressi, 
G. Carugno, R. Onofrio and G. Ruoso, Phys. 
Rev. Lett. 88, 041804 (2002); R. S. Decca, 
D. Lopez, E. Fishbach, G. L. Klimchitskaya, 
D. E. Krause, and V. M. Mostepanenko, Ann. 
Phys. 318, 374 (2005). 

3. T. H. Boyer, Phys. Rev. 174, 1764 (1968). 

4. B. Davies, J. Math. Phys. 13, 1324 (1972). 

5. Annals Phys. 104, 300 (1977); ibid 112, 165 
(1978). 

6. K. A. Milton, L. L. DeRaad and 
J. S. Schwinger, Annals Phys. 115, 388 
(1978). 

7. G. Lambiase, V. V. Nesterenko and 
M. Bordag, J. Math. Phys. (1999), 6254 

hep-th/9812059 . 

8. G. Esposito, A. Y. Kamenshchik and 
K. Kirsten, Phys. Rev. D62, 085027-1 (2000) 
|hep-th/0 006220 . 

9. G. Barton, J. Phys. A 37, 1011 (2004). 

10. E. M. Lifshitz, Zh. Eksp. Teor. Fiz. 29, 94 
(1956); 

11. P. W. Millonni, The Quantum Vacuum (Aca- 
demic, New York, 1993). 



13 



12. K. A. Milton, J . Phys. A 37, R209 (2004) 
|hep-th/0406024| . 

13. R. B. Balian and C. Bloch, Ann. Phys. (NY) 
60, 401 (I970);ibid 63, 592 (I97l);ibid 64, 
271 (1971) Errata in ibid 84,559 (1974); ibid 
69, 76 (1972);ibid 85, 514 (1974). 

14. W. Lukosz, Physica 56, 109 (1971). 

15. J.S. Dowker and R. Critchley, J. Phys. A9, 
535 (1976). 

16. J. Ambjorn and S. Wolfram, Annals Phys. 
147, 1 (1983) 

17. A. Edery, J. Phys. A 39, 685 (2006). 

18. 1. L. DeRaad, Jr. and K. A. Milton, Ann. 
Phys. (N.Y.) 136, 229 (1981). 

19. P. Gosdzinsky and A. Romeo, Phys. Lett. B 
476, 369 (1998) |hep~th/9809l"99| . 

20. I. Brevik and G.H. Nyland, Ann. Phys. 230, 
321 (1994). 

21. K.A. Milton, A.V. Nesterenko and V.V. 
Nesterenko, Phys. Rev. D 59, 105009 (1999). 

22. I. Klich and A. Romeo, Phys. Lett. B 476, 
369 (2000). 

23. V. V. Nesterenko and I. G. Pirozhenko, 
J. Math. Ph ys. 41, 4521 (2000) 

hep-th/9910097 . 

24. R. Balian and B. Duplantier, Recent Devel- 
opments in Gravitational Physics, Institute of 
Physics Conference Series 176, Ed. Ciufiolini 
et al, Oct. 2004 quant-ph/0408124 . 

25. M. Schaden and L. Spruch, Phys. Rev. A 58, 
935 (1998). 

26. M. C. Gutzwiller, Chaos in Classical Quan- 
tum Mechanics (Springer, Berlin, 1990) with 
references to the author's earlier work. 

27. O. Schroder, A. Scardicchio, R. L. Jaffc, Phys. 
Rev. A 72, 012105 (2005). 

28. A. Bulgac, P. Magierski and A. Wirzba, Phys. 
Rev. D73, 025007 (2006). 

29. M.V. Berry and M. Tabor, 
Proc.Roy.Soc. London A349, 101 
(1976); ibid J. Phys. A10, 371 (1977); 
Proc.Roy.Soc.London A356, 375 (1977). 

30. F. Mazzitelli, M.J. Sanchez, N.N. Scoccola, 
J. von Stecher,Phys. Rev. A 67, 013807 
(2003). 

31. M. Schaden, Phys. Rev. A 73,042102 (2006). 

32. C. S chubert, Ph ys. Rept. 355, 73 (2002) 
|hep-th/0101036| ; 



33. H. Gies and K. Langfeld, Int.J.Mod. Phys. 
A 17, 966 (2002) [hep-th/011 2198 ;H. Gies, 
K. Langfeld and L . Moyaerts, J HEP 
0306, 018 (2003) hep-th/0303264 ;ibid, 

hep-th/0311168 ; H. Gies and K. Kling- 
muller, preprint quant-ph/0601094 

34. H. Gies and K. Klingmuller, Proceed- 
ings of QFE XT'05 (Barcelona, Sept. 
'05) 



hep-th /0511092 
quant-ph/0601094. 



ibid, preprint 



35. T. Emig, A. Hanke, R. Golestanian and 
M. Kardar, Phys. Rev. Lett. 87, 260402 
(2001); R. Buscher and T. Emig, Phys. 
Rev. Lett. 94, 133901 (2005); T. Emig, 
R. L. Jaffe, M. Kardar and A. Scardicchio, 
preprint cond-mat/0601055 

36. G. Barton, J. Phys. A 34, 4083 (2001). 

37. I. Cavero-Pelaez and K. A. Milton, 
Ann. Phys. ( N.Y.) 320, 108 (2005) 
|hep-th/0412135| . 

38. A. Romeo and K.A. Milton, Phys. Lett. B 
621, 309 (2005) |hep-th/0504207| . 

39. A. Einstein, Verh. Dtsch. Phys. Ges.19, 82 
(1917). 

40. J. B. Keller, Ann. Phys. (NY)4, 180 (1958); 
ibid in Calculus of Variations and its Appli- 
cations (Am. Math. Soc, Providence, 1958), 
p.27;ibid, J. Opt. Soc. Am. 52, 116 (1962). 

41. V. P. Maslov, Theorie des Perturbations 
et Methodes Asymptotiques (Dunod, Paris, 
1972) and V. P. Maslov and M. V. Feodiriuk, 
Semi- Classical Approximations in Quantum 
Mechanics (Reidel, Boston, 1981). 

42. D. Deutsch and P. Candelas, Phys. Rev. D 
20, 3063 (1979); P. Candelas, Ann. Phys. 
(NY) 143, 241 {1982);ibid,167, 257 (1986). 

43. J.S. Dowker, Ann. Phys. (N.Y.) 62, 361 
(1971). 

44. P. Chang and J.S. Dowker, Nucl. Phys. B395, 
407 (1993). 

45. J.S. Do wker, preprint |hep-th/0510248| ibid, 
preprint math.sp/0601334 . 

46. C. M. Bender, K. A. Milton, Phys. Rev. D 
50, 6547 (1994); K. A. Milton, Phys. Rev. D 
55, 4940 (1997). 

47. A. A. Sahari an, Preprint IC/2000/14, 
|hep-th/ 0002239 

48. F. Bernasconi, G.M. Graf and D. Hasler, 



11 



Ann. Inst. H. Poincare 4, 1001 (2003) 
|math-p h/0302035|. 

49. M. Bordag,K. Kirsten, D. Vassilevich, Phys. 
Rev. D59, 085011 (1999). 

50. N. Graham, R. L. Jaffe, V. Kle- 
mani, M. Quandt, M. Scandurra and 
H. Weigel, Nucl. Phys. B645, 49 (2002); 
R. L. Jaffe, MP Conf. Proc. 687, 3 (2003) 



arXiv:hep-th/0307014 . 



51. K. A. Milton, Phys. Rev. D68,065020 (2003); 
ibid J. Phys. A37, R209 (2004). 

52. A. A. Saharian and A. S. Taroyan, 
|hep-th/0603144| 

53. A. Scardicchio, R. L. Jaffe, Nucl. Phys. 704. 
552 (2005);i6zd, preprint quant-ph/0507042. 

54. J.B. Keller, J. Opt. Soc. Amer. 52, 116 
(1962); J.B. Keller and R.M. Lewis, Asymp- 
totic Theory of Wave Propagation and 
Diffraction, unpublished book manuscript in 
2 Volumes, (Courant Institute, New York 
University, 1970);J.B. Keller, J. Appl. Phys. 
30, 1452 (1952);B.R. Levy and J. Keller, 
Commun. Pure Appl. Math. 12, 159 (1959); 
H.M. Nussenzveig, Ann. Phys. 34, 23 
(1965);H. Primack, H. Schranz, U. Smilansky 
and I. Ussishkin, J. Phys. A:Math. Gen. 30, 
6693 (1997). 

55. M. Schaden and L. Spruch, Ann. Phys. (N.Y.) 
313, 37 (2004). 

56. G. Vattay, A. Wirzba, P.E. Rosenqvist, Phys. 
Rev. Lett. 73 (1994) 2304. 

57. A. Romeo, private communication to 
K. A. Milton (reproduced in App. B of|2l]L 

58. M. Bordag and I. G. Pirozhcnko, Phys. Rev. 
D 64, 025019 (2001). 

59. V.V. Nesterenko and I. G. Pirozhenko, Phys. 
Rev. D 60, 125007 (1999). 

60. E. A. Power, Introductory Quantum Electro- 
dynamics (Elsevier, New York, 1964), Ap- 
pendix I. 

61. L. Brown and G. J. Maclay, Phys. Rev. 184, 
1272 (1969); M. Schaden and L. Spruch, 
Phys. Rev. A 65, 034101 (2002). 

62. J. S. Hoye, I. Brevik, J. B. Aarseth, K. A. Mil- 
ton, quant-ph/0506025 

63. P.M. Neto, A. Lambrecht and S. Reynaud, 
|quant-ph/05050"86| 



